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The experimental concept of a search for a long-range coupling between rubidium (Kb) nuclear 
spins and the mass of the Earth is described. The experiment is based on simultaneous measurement 
of the spin precession frequencies for overlapping ensembles of *^Rb and *^Rb atoms contained 
within an evacuated, antirelaxation-coated vapor cell. Rubidium atoms are spin-polarized in the 
presence of an applied magnetic field by synchronous optical pumping with circularly polarized 
laser light. Spin precession is probed by measuring optical rotation of far-ofi'-resonant, linearly 
polarized laser light. Simultaneous measurement of *^Rb and *^Rb spin precession frequencies 
enables suppression of magnetic-field-related systematic effects. The nuclear structure of the Rb 
isotopes makes the experiment particularly sensitive to anomalous spin-dependent interactions of 
the proton. Experimental sensitivity and a variety of systematic effects are discussed, and initial 
data are presented. 
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I. INTRODUCTION 



The connection between quantum theory and general 
relativity is one of the most important unsolved mys- 
teries of modern physics, a mystery exacerbated by the 
dearth of experiments probing the rare intersections be- 
tween these two theories. One intersection between quan- 
tum effects and gravity that does offer potential for ex- 
perimental tests is the question of how intrinsic spins 
interact with gravitational fields (see, for example, the 
review 'l| ) . According to general relativity, a purely ten- 
sor theory, the intrinsic spin of a particle is unaffected by 
the local gravitational field 0tS|- However, in extensions 
of general relativity based on a Riemann-Cartan space- 
time instead of a Riemann geometry, the gravitational 
interaction is described by a torsion tensor which can 
generate heretofore undetected spin-mass and spin-spin 
interactions |7l-[l0|. 

In terms of quantum field theory phenomenology, 
the torsion tensor describes new scalar-pseudoscalar 
and vector-pseudovector gravitational interactions, cor- 
responding, respectively, to spin-0 and spin-1 gravitons 
in addition to the usual spin-2 graviton associated with 
the tensor nature of standard gravity |lll - [l3| . New spin-1 
and spin-0 partners of the usual spin-2 graviton also nat- 
urally arise in theoretical attempts to unify gravity and 
quantum mechanics, such as string theory and M-th eory 
p^l - llq , especially in the context of supersymmetry [l^l . 
It has recently been noted that a massless or nearly mass- 
less spin-0 component of gravity manifests as dark energy 
over cosmological distances |18l - [21| . The pseudoscalar 
component of such a field leads to an interaction that 



has the nonrelativistic form 12 ill: 



Hg = k-cr g 



(1) 
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where fc is a dimensionless parameter setting the scale 
of the new interaction, h is Planck's constant, cr is the 
intrinsic spin of the particle in units of h, g is the Earth's 
gravitational field, and c is the speed of light. If the 
strength of the pseudoscalar coupling is the same as that 
of the tensor component of gravity. A: w 1 [21]. 

The Hamiltonian Hg in Eq. ([1]) manifestly violates 
the equivalence principle for intrinsic spins, offering a 
mechanism by which a gravitational field can be distin- 
guished from an accelerating reference frame. Further- 
more, if Hg y^ 0, gravity would violate parity (P) and 
time-reversal (T) symmetries, and the spin-gravity cou- 
pling would be a source of additional CP-violation that 
might explain the observed matter-antimatter asymme- 
try of the universe [22| ■ The possibility that gravity may 
violate discrete symmetries has led many authors over 
the past fifty years [3, [23l - [30| to consider the possibility 
of such an interaction, which would imply the existence of 
a gravitational dipole moment (GDM) khcr/c for elemen- 
tary particles. One could envision a GDM as a separa- 
tion between the center of inertial mass and the center of 
gravitational mass by a distance kh/mc, where m is the 
particle mass. If the dimensionless coupling constant k 
is of order unity, then the separation between the centers 
of inertial mass and gravitational mass is on the order of 
the Compton wavelength. Generally, theoretical models 
of gravity that accommodate such an interaction predict 
that fc ^ 1 [3, I23I430II. However, to date, the most sen- 
sitive searches [Slj, Wl\ for such a spin-gravity coupling 
have set limits fc < 10, still an order of magnitude away 
from the most theoretically interesting region of param- 
eter space. The central goal of the experiment described 
in the present paper is to probe spin-gravity interactions 
of the type described by Eq. ([T]) at the fc ~ 1 level. 

There are several experimental consequences of the ex- 
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FIG. 1: Existing experimental constraints (at the 2-a level) on nucleon monopole-dipole (scalar-pseudoscalar) couplings 
IgpQsl /he as a function of the range A of the interaction {gp and gs are the pseudoscalar and scalar coupling constants, 
respectively). Direct experimental constraints for the neutron are from Youdin et al. (1996) [33|] at the laboratory-scale range 
and from Venema et al. (1992) [3l|] for the earth-scale range (excluded parameter space shaded blue); constraints for the 
proton are from the experiment of Wineland et al. (1991) [34| (excluded parameter space shaded purple). Astrophysical con- 
straints for baryon couplings (excluded parameter space shaded green) are from the recent analysis of Raffelt (2012) 35]. The 
nominal coupling strength for a scalar-tensor theory of gravity corresponding to fc ~ 1 is represented by the blue line at the 
Igpflsl /he « 10~^® level. Constraints on monopole-dipole couplings of the electron spin to the mass of the earth (obtained by 
the University of Washington torsion pendulum experiment [32l, not shown) are similar to the constraints on the neutron spin 
from Ref. [3l[. The potential sensitivity (dashed red line) of our proposed search could improve upon existing experimental 
limits on long-range (A > 10^ m) monopole-dipole couplings in general by an order of magnitude and for the proton spin in 
particular by three orders of magnitude. 



istence of a GDM for an elementary particle. The inter- 
action described by Eq. ([T]) leads to a gravity-induced 
splitting AE of the energy levels for spins oriented par- 
allel and anti-parallel to g: 



AE^2k^ 



4fc X 10^23 gy _ 



(2) 



In addition, an interaction such as that described by 
Eq. dl]) generates a torque on spins immersed in a grav- 
itational field, leading to spin precession about the axis 
of the local gravitational field with a frequency 



n„ 



kg 



-^ = -1 w fc X IQ-'^ Hz 

ZTT TTC 



(3) 



This spin precession not associated vifith the magnetic 
moments of the particles is the signature of the P- and 
T-violating spin-gravity coupling that is being searched 
for in our present experiment. 

In principle, the effect of a spin-gravity coupling is in- 
distinguishable from the interactions resulting from any 



heretofore undiscovered force-mediating pseudoscalar or 
vector particle. In the literature, such new interactions 
are commonly parameterized using the Moody- Wilczek 
formalism for spin-0 particles [36|, recently extende d by 
Dobrescu and Mocioiu to include spin-1 particles |37| . 
In this context, the spin-gravity (or spin- mass) inter- 
action searched for in our proposed experiment can be 
interpreted as a monopole-dipole coupling. In general, 
the monopole-dipole coupling strength can be different 
for different elementary particles — our proposed experi- 
ment is primarily sensitive to monopole-dipole couplings 
of the proton, whereas the previous best jsxperimental 
limits are for neutron [31| and electron (33] couplings. 
Consequently, our experiment has the potential to im- 
prove experimental constraints on monopole-dipole cou- 
plings of the proton by orders of magnitude compared to 
the best previous limit |34| . 

Figure [1] presents a parameter exclusion plot showing 
existing direct experimental limits on monopole-dipole 
interactions of nucleons at various length scales, as well 
as astrophysical constraints inferred from the duration 



of the supernova SN 1987A neutrino burst in combina- 
tion with data from searches for anomalous monopole- 
monopole forces [3a]- (It should be noted that there are 
significant uncertainties related to dense nuclear matter 
effects in the analysis of the SN 1987A neutrino burst 
[381 Isgj.l The blue horizontal line at the bottom of the 
plot shows the strength of the monopole-dipole interac- 
tion corresponding to fc w 1 . The dashed red curve shows 
the projected sensitivity of our experiment. Note that as- 
trophysical limits for electron couplings, which are based 
on star cooling, are three orders of magnitude more re- 
strictive than astrophysical limits on nucleon couplings. 
Thus the astrophysical limits for electrons just reach the 
/c « 1 regime, and are comparable to the projected sen- 
sitivity of our experiment. Of interest are observations 
that the white-dwarf luminosity function fits better with 
a small amount of anomalous energy loss at a level that 
would correspond to a pseudoscalar interaction with cou- 
pling strength corresponding to /c « 1 [40| . Furthermore, 
the observed period decrease of the pulsating white dwarf 
G117-B15A also favors some amount of extra cooling [41j . 
In the present work, we describe a dual-isotope ru- 
bidium (Rb) comagnetometer well-suited for searching 
for a long-range monopole-dipole coupling between pro- 
ton spins and the mass of the earth. We show that the 
dual-isotope Rb comagnetometer as designed can achieve 
sufficient statistical sensitivity and rejection of known 
systematic errors so that present experimental limits on 
anomalous monopole-dipole couplings of the proton spin 
can be improved by orders of magnitude. The basic con- 
cept of our experiment is to use synchronous laser op- 
tical pumping to generate transverse spin polarization 
of Rb atoms contained in an antirelaxation-coated cell 
[42:j , and then employ off-resonant laser light to simulta- 
neously measure the spin precession frequencies of *^Rb 
and ^''Rb atoms in the presence of a magnetic field B. 
The ratio of the difference between the Rb precession 
frequencies divided by their sum. 
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will be measured for a range of different magnetic fields. 
Measurement of the ratio 3? eliminates or reduces sev- 
eral common-mode sources of noise and systematic error. 
Taking the difference A3? between 31 for B parallel with 
g and anti-parallel with g yields a signal proportional 
to the spin precession frequency caused by nonmagnetic 
interactions. In this configuration, the valence electron 
spin of the Rb atoms effectively serves as an accurate 
comagnetometer for the Rb nuclear spins. 

Our early work on this experiment is described in 
Ref. [43]. In this early version of the experiment, a sin- 
gle linearly polarized laser beam, frequency-modulated at 
an integer multiple of the precession frequency, was used 
to measure nonlinear magneto-optical rotation (NMOR) 
induced by the spin precession of the Rb atoms [44| - |47| . 
While the NMOR measurements achieved our target sta- 
tistical sensitivity to the Rb spin precession frequencies. 



a subtle systematic effect involving light shifts |43| . |48| re- 
quired us to modify our experimental approach. The ef- 
fect, known as alignment-to-orientation conversion, arises 
from the combined action of the magnetic field and op- 
tical electric field. Alignment-to-orientation conversion 
evolves spin polarization aligned along the linear polar- 
ization axis into spin polarization oriented along the light 
propagation (magnetic field) direction. Atomic spins ori- 
ented along the light propagation direction generated el- 
lipticity of the light field which in turn produced unac- 
ceptably large light-power- and magnetic-field-dependent 
shifts of the spin precession frequencies. Our present ex- 
periment circumvents these systematic effects by tempo- 
rally separating pump and probe stages, using unmod- 
ulated probe light, and by a choice of experimental ge- 
ometry where the probe beam propagates in a direction 
orthogonal to the magnetic field (Sec. W} . 



II. EXPERIMENTAL CONCEPT 

As noted in Sec. U a long-range spin-mass or spin- 
gravity interaction can be parameterized in terms of a 
GDM K via /t = x^^i where x = kh/c is the "gyro- 
gravitational ratio" for the particle. The x for Rb atoms 
can be calculated in terms of Xe and Xp-i the gyro- 
gravitational ratios for the electron and proton, respec- 
tively, using the shell model to describe the nuclei. In 
a given ground-state hyperfine level with total angular 
momentum F, the atomic GDM Katom(-F') is 

/«atom(F) = X.ton.{F)F (5) 

^S^e^Fl {I-F) 

F{F + 1) ^ ' F{F + 1) ^ ' ' 

where Se is the electron spin, I is the nuclear spin, and 
Xnuci is the nuclear gyro-gravitational ratio. In the nu- 
clear shell model [43] , both Rb isotopes have valence pro- 
tons, and so the nuclear spin, magnetic moment, and nu- 
clear GDM are, to a good approximation, due entirely to 
the proton. The nuclear GDM Knuci is thus given by 



Xnnc\I ■ 



(S„ ■ I) 



{I)Xp , 



(6) 



/(/+1) 

where Sp is the proton spin and we have assumed, as 
do most theoretical models [3, [2l|, |23| - |30| , that there is 
no contribution from orbital angular momentum. Rele- 
vant parameters based on the above considerations are 
presented for the two Rb isotopes in the probed ground- 
state hyperfine levels in Table HI 

Ignoring temporarily other causes of spin precession, 
the spin-precession frequencies for *^Rb and ^^Rb in the 
presence of the magnetic field B and Earth's gravita- 
tional field g are 
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(8) 



TABLE I: Parameters determining gyro-gravitational ratios x ^'iid gyromagnetic ratios 7 = gFMo for *^Rb and *^Rb in the 
ground-state hyperfine levels of interest. 



Atom Ground electronic state Nuclear spin Total angular momentum g-factor Proton state Xatom 
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7^ = 2 


9F = 1/2 
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where 4> is the angle between B and g, 735 and 737 are the 
gyromagnetic ratios (7 = gpfJ-o, where gp is the Lande 
5-factor and /zo is the Bohr magneton) and the light prop- 
agation direction is along B. In the above we neglect con- 
tributions to the spin-precession frequency second-order 
in g. To analyze the data, we construct the following 
ratio: 
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To first order assuming jB ^ Xeg, Xp9 a-^id neglecting 
the elTects of the nuclear magnetic moments, we have 
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f 787 - 785 
V787 + 785 
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Xp9 cos ' 
fJ-oB 



(10) 



where 3?+ is for positive B and !K_ is for negative B 
(relative to g). There is first-order cancelation of the 
effects of an electron GDM in the ratio 5i, and near unity 
sensitivity to the effects of a proton GDM. Measuring 
AJi = D?+ — !K_ yields a signal proportional only to the 
proton GDM: 



AJIk4.12\ 



787 — 785 
'V 787 +785 



Xp9 cos t 
fioB 



(11) 



The first-order cancelation of the electron GDM contri- 
bution to A'R is a result of the fact that we measure spin 
precession in the 7^ = / -I- 1/2 ground state hyperfine 
level for both isotopes, so electron couplings contribute 
in nearly identical ways to the measured values of 1^85 
and $787- 



III. EXPERIMENTAL SETUP 

A schematic diagram of the experimental setup used 
to carry out simultaneous measurement of ilss and r287 
is shown in Fig. [2] At the heart of the experiment is 
a natural isotopic mixture of Rb vapor (72.2% ^^Rb, 



27.8 



87 



Rb) contained within an evacuated (residual 



pressure « 10~^ torr) spherical alkene-coated glass cell 
(diameter = 5 cm). The alkene coating is 1-nonadecene 
[CH2 — CH(CH2)i6 — CH2] and the cell wa s prepared ac- 



cording to procedures outlined in Ref. 42|. The partic- 
ular cell we are using was measured to have longitudi- 
nal spin relaxation times Ti « 5 s limited by exchange 



of atoms between the spherical bulb of the cell and the 
stem which contains the Rb reservoir. Under typical op- 
erating conditions, the spin relaxation rate due to wall 
collisions is significantly smaller than the relaxation rate 
due to spin-exchange collisions between the Rb atoms. 

The vapor cell is mounted inside a frame manufactured 
of HDPE (High Density Polyethylene) which is fit inside 
the innermost layer of a five-layer magnetic shield (man- 
ufactured by Amuneal Inc.) made of a 1-mm thick high- 
permeability alloy, annealed in a hydrogen atmosphere. 
Each layer of the shield consists of a cylindrical center 
piece and two removable end caps. The outer layers of 
the shield are spaced by styrofoam (polymerized in place) 
and the innermost layer is spaced by melamine foam to 
reduce acoustic noise. Four ports for access to the in- 
side of the shields are available on the cylindrical pieces 
and one port is available on each end cap. The shield- 
ing factor of the entire five-layer magnetic shield system 
was measured for a nearly identical design to be better 
than lO'' [5i|. The foam spacing between the shield lay- 
ers provides thermal insulation in addition to mechanical 
support. The temperature of the innermost shield layer 
is stabihzed at 30°C by a J-KEM Model 210 temperature 
controller using a T-type thermocouple attached to the 
inner surface of the shield layer for temperature measure- 
ment and resistive heating with a twisted pair of wires 
wrapped about the outside of the innermost shield layer. 
Stabilizing the shield temperature serves two functions: 
(1) it reduces temperature- related drifts of residual mag- 
netic fields from the innermost shield and (2) it provides 
a stable, elevated temperature environment for the Rb 
cell yielding vapor densities of w 2 x 10^° atoms/cm^. 

A system of nine separate coils are wound in grooves 
cut into the frame mounted inside the innermost layer of 
the shield. The system of coils was designed to provide, 
over the volume of the Rb vapor cell, uniform magnetic 
fields in three orthogonal directions {B^, By, and B^), 
linear magnetic field gradients in five directions (dBx/dx, 
dBz/dz, dBx/dz, dBy/dz, dBy/dx), and a quadratic gra- 
dient along the shield axis (d^Bz/dz^). As a consequence 
of Maxwell's equations (V • B = and V x B = 0), con- 
trol over the five linear magnetic field gradients is suf- 
ficient to provide compensation of all nine possible lin- 
ear gradients. Based on computer modeling (using the 
Amperes program from Integrated Engineering Software 
Inc.), the uniformity of the magnetic fields and linear- 
ity/quadracity of the field gradients generated by the coil 
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FIG. 2: Right-hand side: schematic diagram of the experimental setup used to measure the spin precession frequencies of *^Rb 
and ^'^Rb and to search for a long-range spin-mass (spin-gravity) coupling (P = linear polarizer, M = mirror, PBS — polarizing 
beamsplitter, A/4 = quarter wave plate). Picture, upper left: schematic of the apparatus used to align the magnetic shield 
system and magnetic field B along the Earth's rotation axis VLe, and geometrical relationship to the local gravitational field g 
(this geometry is chosen to control systematic errors related to the earth's rotation, see Sec. lVB]l . 



system is at a part per thousand over the cell volume for 
typical applied currents. It should be noted that effects of 
uncompensated magnetic-field grad ients are significantly 
reduced by motional averaging [50| (effects are quadratic 
in the the magnitude of the gradient). The coils are in 
series with a set of ultra-stable, low temperature coeffi- 
cient (low TC) resistors (Caddock Type USF 200 Series, 
zero nominal TC with TC < 2 ppm/K). The voltage for 
the Bz coil is supplied by a precision DC voltage source 
(Krohn-Hite Model 523 calibrator, stability ±1 ppm) and 
voltages for the coils controlling Br^, By, and field gra- 
dients are computer generated with a digital-to-analog- 
converter (DAC, National Instruments PCI-6733). 

In order to measure rjgs and fisTi a- system of shut- 
ters (Stanford Research SR474) is used to implement a 
temporally separated pump/probe measurement scheme. 
During the optical pumping stage (duration = 1 s), Rb 
atoms are illuminated by two collinear, circularly po- 
larized pump beams propagating along — y (orthogonal 
to B which is along z), one tuned to the center of the 
Doppler-broadened ®^Rb D2 F = 3 — )■ i^' resonance and 
the other tuned to the center of the Doppler-broadened 
s^Rb Dl F = 2 ^- F' = 1 resonance (F, F' are the to- 
tal atomic angular momenta of the ground and excited 



states, respectively). The 780- nm D2 pump beam is 
generated by a distributed feedback laser diode (EYP- 
DFB-0780-00080-1500-TOC03 from Eagleyard Photon- 
ics) and the 795-nm Dl pump beam is produced by 
a tunable external-cavity diode laser (Toptica DLIOO). 
The pump beams are amplitude-modulated at frequen- 
cies close to the respective Larmor frequencies of the 
isotopes using electro-optic modulators (EOMs, Thor- 
LABs EO-AM-NR-Cl, not shown in Fig. [2]) placed be- 
tween crossed calcite linear polarizers. The duty cycle 
for both pump beams is 20%; during the period when 
the EOMs transmit the pump light, the power of the 
D2 pump beam incident on the Rb atoms is k, 55 /^W 
and the power of the Dl pump beam is « 150 /zW. 
These parameters were chosen to maximize the trans- 
verse spin polarization for both Rb isotopes. This syn- 
chronous optical pumping generates atomic spin polar- 
ization transverse to B in both isotopes precessing at 
their respective Larmor frequencies. A third, collinear, 
circularly polarized re-pump beam tuned to the center 
of the Doppler-broadened ^^Rb D2 F = 1 — > F' reso- 
nance transfers *^Rb atoms pumped into the unobserved 
F = 1 ground-state hyperfine level back into the F = 2 
hyperfine level (which, taking into account natural iso- 



topic abundances, yields approximately equal signals for 
both isotopes). The 780- nm re-pump beam is produced 
by tunable external-cavity diode laser (New Focus TLM 
7000), with power « 750 fiW. The diameters of the pump 
and re-pump beams are « 2 mm. 

During the optical probing stage (duration = 1 s), a 
shutter blocks the pump and re-pump beams, and shut- 
ters open to allow a linearly polarized probe laser beam 
to propagate along — y through the vapor and into a 
polarimeter. The 780-nm D2 probe beam is produced 
by another tunable external-cavity diode laser (Toptica 
DLIOO). The ^^Rb and ^^Rb precession frequencies, figs 
and rig?, are measured by observing optical rotation 
of the probe light. The frequency of the probe beam 
is tuned w 3 GHz below the center frequency of the 
Doppler-broadened ^^Rb D2 F = 2 —>■ F' resonance, the 
power is ~ 200 /iW, and the beam diameter is « 2 mm. 
Prior to entering the vapor cell, the probe beam passes 
through an antireflection-coated Glan Thomson linear 
polarizer (calcite, extinction ratio 5 x 10^ : 1). After ex- 
iting the vapor cell, the beam is analyzed by a polarime- 
ter consisting of a WoUaston prism polarizing beamsplit- 
ter (calcite, extinction ratio 10^ : 1) whose output rays 
are detected with a balanced photoreceiver (New Focus 
Model 2307). The signal from the photoreceiver is sent to 
a preamplifier (Stanford Research Systems SR560) and 
then recorded on computer using an analog-to-digital 
converter (National Instruments PCIe-6361) using a rou- 
tine written in Lab VIEW. The time base for the data 
acquisition is provided by a 10 MHz signal from a Rb 
atomic frequency standard (Stanford Research Systems 
SIM940, short-term stability < 2 x lO'^^ in 100 s) that is 
GPS-disciplined with a 1 PPS signal (from a Communi- 
cation Navigation Surveillance Inc. CNS Clock II, with 
a long-term accuracy better than a part in « 10^^). The 
accurate time base ensures that flg^ and figy can be mea- 
sured at the 10~® Hz level over a long period of time for 
data averaging. The pump and re-pump lasers are fre- 
quency stabilized using dichroic atomic vapor laser locks 
(DAVLLs) [53,l53|. The probe, pump, and re-pump beam 
spectral purities are monitored with Fabry-Perot inter- 
ferometers (ThorLABs SA200-5B) and the light powers 
of each beam transmitted through separate uncoated Rb 
reference cells (natural isotopic mixture) are monitored 
to ensure that the lasers remain properly tuned. (The 
laser frequency locking and diagnostics setups are not 
shown in Fig. [21) 

The picture in the upper left corner of the experimental 
setup diagram (Fig. [2]) depicts the magnetic shield mount 
used for mechanical alignment of the shield axis z along 
the Earth's rotation axis CIe , which is important for con- 
trol of a systematic error related to the Earth's rotation 
(Sec. IV B|) . The outermost shield layer is held in place 
with an aluminum frame, which is bolted to precision tilt 
and rotation stages (Newport TGN160 and UTR120, re- 
spectively) attached to an optical breadboard tilted from 
horizontal by an angle approximately equal to the lati- 
tude of the laboratory (37° 39' 24" N). By surveying the 



laboratory (using Google Earth as well as GPS signals) 
and using an alignment laser propagating along the shield 
axis (z) with a path length of « 3 m, we are able to me- 
chanically align the shield axis with the Earth's rotation 
axis to within 0.3°. Prior to measurement of Qgc, and 
rig7, B is carefully aligned along the z-axis using a laser 
beam split off from the 780-nm probe beam by measuring 
NMOR with frequency modulated light 43, 44, 54, 55,] . 
The accuracy of the alignment of B along the light propa- 
gation direction k using NMOR ^] is much greater than 
the mechanical alignment accuracy of the shield axis, so 
the alignment of B parallel with (^e is achieved with an 
uncertainty of « 0.3°. 



IV. INITIAL DATA AND PROJECTED 
STATISTICAL SENSITIVITY 



Sample data acquired during the probe sequence are 
shown in Fig. [3l The upper plot shows the optical ro- 
tation signal acquired in the time domain for an applied 
field of B sa 7.1443 mG oriented in the z direction, and 
the lower plot shows the absolute value of the Fourier 
transform of the data set (carried out using a data anal- 
ysis routine written in Mathematica) . In the frequency 
domain, distinct resonant peaks in the Fourier transform 
can be identified and correspond to fls5 and figy. The 
dominant contribution to these signals is from atoms in 
the resonantly pumped ^^Rb F = 3 ground state and the 
^^Rb F = 2 ground state, since the pump and re-pump 
laser beam parameters are optimized for transverse spin 
polarization of these states and the probe beam is tuned 
closest to optical resonance with transitions from these 
states. 

In the time domain a slow beating is observed in the 
optical rotation signal. The additional frequency com- 
ponent responsible for the beating arises from spin pre- 
cession of ^^Rb atoms in the F = 2 ground state, which 
has a slightly different Lande ^-factor magnitude than 
the F = 3 ground state because of the nuclear magnetic 
moment (this effect is discussed in detail in Sec. IV C|) . 
The ^^Rb F = 2 ground state is slightly polarized by 
nearly synchronous re-population pumping: a fraction of 
the ®^Rb atoms optically pumped from the F = 3 ground 
state spontaneously decay from the 5 ■^P3/2 excited state 
back to the F = 2 ground state; the absolute values of 
the Lande g-factors are close enough in value that for 
sufficiently small magnetic fields a detectable transverse 
spin polarization in the ®^Rb F — 2 ground state can be 
created during the pumping stage. (^^Rb atoms in the 
-F = 1 ground state contribute a much smaller amplitude 
signal for a variety of reasons, discussed in Sec. IV CI ) 

Data analysis is carried out by fitting subsets of 
the Fourier transformed optical rotation signal centered 
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FIG. 3: Upper plot: sample time-domain data from 
pump/probe measurement of Rb spin precession using opti- 
cal rotation. Inset in upper right corner shows a subset of the 
data of total duration 10 ms (highlighted by the red box on 
the complete data set). In the inset, fast beating is observed 
between sinusoidal signals with frequencies figs and ilg?. The 
applied magnetic field corresponds to |B| ~ 7.1443 mG. Slow 
beating can be observed in the main time-domain plot be- 
tween a dominant signal from *^Rb atoms in the F = 3 ground 
state hyperfine level and a smaller-amplitude signal from *^Rb 
atoms in the F — 2 ground state hyperfine level (the absolute 
value of the Lande g-factors of the two ground state hyper- 
fine levels differ due to the nuclear magnetic moments, see 
Lower plot: Fourier transform of the time-domain 



sec. rvq) 

data, showing resonances at figs (lower-frequency peak) and 
fls7 (higher-frequency peak) . 



around fJgs ^^nd Jlgy to a Lorentzian function S{io): 



Siu:) 



a 



\U+{^)\ 



I3{uj -n)/r 



I 



(12) 



where a and /3 are the amphtudes of the imaginary and 
real components of the signal, respectively, uj is the fre- 
quency, n is the resonant spin-precession frequency, and 
F is the resonance width (corresponding to the spin re- 
laxation rate). To account for spin precession of ^^Rb 
atoms in the F — 2 ground state, additional real and 
imaginary Lorentzian components of the signal can be 
included in the fitting function for fls5- Under typical 
operating conditions, F/(27r) w 1 Hz. 



Fits to the data demonstrate a statistical sensitivity to 
the spin precession frequency f2 of 



6n 

2^ 



100 fiTiz 



(13) 



for a 1 s measurement. Ultimately, the shot-noise- limited 
(SNL) sensitivity (Jf^sNL of a spin-polarized atomic sam- 
ple to precession frequencies is determined by the to- 
tal number of atoms TV and the relaxation rate Fj-oi 
of the atomic spin polarization (for measurement times 
r > Frei"' Mm)- 



sa 



SNL 



Trol 

'Nt 



(14) 



Under our experimental conditions, TV sa 10^^ atoms and 
Fi-oi ~ 27r X 1 Hz (limited by spin-exchange collisions), 
yielding (JfJsNL ~ 27r x 0.4 ^Hz for t = 1 s. This suggests 
that the measurement is not presently shot-noise-limited, 
and further reduction in technical noise would permit 
even better statistical sensitivity. 

We are in the process of making several modifications 
to the apparatus in order to improve the sensitivity. 
Noise from the balanced photoreceiver presently exceeds 
the photon-shot-noise limit by a factor of « 5, and there- 
fore we are upgrading the balanced photoreceiver in order 
to achieve photon-shot-noise-limited polarimetry. Opti- 
mization of the probe light power and detuning should 
enable further improvement in sensitivity, since the pho- 
ton shot-noise limit presently exceeds the atomic shot- 
noise limit by over an order of magnitude. 

Although comagnetometry significantly reduces 
magnetic-field-related noise and systematic effects 
from acquisition-to-acquisition, in our configuration 
it does not improve the statistical uncertainty for a 
single acquisition. Therefore magnetic field noise can 
degrade the sensitivity. In particular, magnetic field 
noise due to thermal Johnson currents in the innermost 
mu-metal shield is estimated to contribute noise at the 
100 /iHz/vHz level |6ll - l63l |. We are in the process of 
replacing the innermost mu-metal shield with a non- 
conducting ferrite shield, which has been demonstrated 
to reduce thermal magnetic field noise [64|. 

We expect that these improvements to our apparatus 
should enable sensitivity to atomic spin precession at the 
10 /iHz level for a 1 s measurement. Collecting data 
for sa 10^ s would then yield a statistical sensitivity of 
w 10~* Hz to anomalous spin-precession, sufficient to 
search for a GDM signal corresponding to fc ^ 1 (Eq. [3]). 



V. SYSTEMATIC EFFECTS 

A. General considerations 

For a general consideration of systematic errors, it is 
helpful to characterize the ways in which additional con- 
tributions to the spin-precession frequencies for ^^Rb and 
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*^Rb beyond Larmor precession can enter the expressions 
for figs and ilg?: and, crucially, the comagnetometer sig- 
nal A3i from which we will extract the GDM coupling. 

One useful way to characterize systematic errors is to 
separate contributions to the spin precession frequencies 
into those that reverse sign when the direction of B is 
changed relative to g (_B-odd terms, fto) and those that 
do not reverse sign (_B-even terms, fig). Depending on 
the orientation of B relative to g, for each isotope we 
obtain two different precession frequencies 



n± =QL+^e± ^o 



(15) 



where ^II is the appropriate Larmor frequency. Some- 
what counter-intuitively, any effect that causes spin pre- 
cession in a fixed sense contributes a i?-odd term. This 
is because reversal of B reverses the sense of Larmor pre- 
cession, and it is the absolute value of the spin preces- 
sion frequency that is measured in the experiment. Thus 
spin-prcccssion due to a GDM coupling is a B-odd term 
contributing to flo, and consequently i3-odd systematic 
errors are not suppressed in the comagnetometer signal 
AIR. Assuming J7e = and a B-odd systematic effect 
adding to the *^Rb precession frequency, we have 



AD? « 9.6 



787 - 785 \ ^o 



787 + 785 / MO-B 



(16) 



Therefore B-odd systematic effects must be suppressed 
or accounted for by other means. 

In the case of perfect magnetic-field reversal, there is 
no contribution of B-even terms to AJl. However, if 
magnetic- field reversal is imperfect by an amount SB, 
B-even systematic effects can lead to a nonzero AJl. As- 
suming flo = and a i3-even systematic effect adding to 
the *^Rb precession frequency, we have 



AJi'. 



g/ 787-785 \ 
V787+785/ 



f7e SB 



(17) 



Thus B-even systematic effects are suppressed relative 
to i?-odd effects by a factor ~ SB/B which in our ex- 
periment can be made < 10^^ by taking advantage of 
the high sensitivity of the setup to magnetic fields (un- 
der typical experimental conditions, the applied magnetic 
field is ~ 10 mG and the magnetometric sensitivity of the 
apparatus is SB ^ 10"^"'^ G/-\/Hz). 

Because our experiment employs a scalar measurement 
scheme (see, e.g., Refs. [53,|60|) where the dominant con- 
tribution to the spin precession frequency is from Larmor 
precession induced by the magnetic field B, we can also 
characterize spin precession in terms of contributions 51 1| 
that add linearly to VIl and contributions n± that add 
in quadrature to Larmor precession: 



where we have assumed that ri||, J7j_ <^ il^. 

Certain systematic errors (such as light shifts) are sup- 
pressed by arranging the experimental geometry so that 
they contribute to fl primarily as D,j_. Of course, imper- 
fections in alignment inevitably mean that there is some 
contribution of such systematic errors to both flu and 



ri|| = f^orr simp W ipQcii , 



n±_ = Ocir COS if 



V 



(20) 
(21) 



1 „ I i 'orr J 

n^^L+ V^.rr + %^ , (22) 

where ilorr is the amplitude of the systematic error and </? 
is the misalignment angle from perfect orthogonality to 
the leading field contribution. Experimentally, mechan- 
ical alignment of the system can in most cases achieve 
at best (y5 ^ 5 X 10^^ rad (0.3°). However, in the case 
of the alignment of B parallel with or orthogonal to the 
light propagation direction (represented by the wave vec- 
tor k) , much better results can be achieved by employing 
nonlinear magneto-optical effects that depend on the an- 
gle between B and k [56[: ip < 10~^ rad can be achieved 
under typical operating conditions. 

Other systematic errors cannot be sufficiently sup- 
pressed using the above experimental geometry (for ex- 
ample, the gyroscopic error introduced by rotation of the 
Earth |3l|; 123: l53] ) ■ In these cases, the experiment is ar- 
ranged so that the error contributes to 57 primarily as 



n 



II- 



ri|| = riorr COS (/3 W 1 
Q,±_ — 57crr sin Lp W iy9ilc 



f 



n. 



,2o 2 

err 



n^n, + [i-^y... + ^ 



(23) 
(24) 
(25) 



While this geometry offers no suppression fiorr, the sys- 
tematic uncertainty in the value of J7orr due to apparatus 
misalignment is quadratically suppressed. If ilon- is in- 
dependently measured with sufficient accuracy, it can be 
subtracted from the data. 

The various systematic effects considered in this sec- 
tion are summarized in Table UTJ which lists their esti- 
mated contribution to A3i. The most significant esti- 
mated source of systematic uncertainty in our experiment 
is the effect of light shifts due to the residual ellipticity 
of the nominally linearly polarized probe beam. 



B. Gyro-compass effect 



Q 



{nL + n^y + in^y 



^r. + ^u 



<^2 



2n, 



(18) 
(19) 



Because the experimental apparatus is attached to the 
Earth, while the atomic spins are decoupled from Earth's 
rotation, the experimental signal is sensitive to the ro- 
tation rate of the Earth, J7E/(27r) w 11.6 /iHz. This 



TABLE II; Estimated contribution of various systematic er- 
rors to AJi for |B| = 7.1443 mG. The atomic shot-noise- 
limited sensitivity of the setup under our experimental con- 
ditions, A'' ~ 10^^ atoms and F^gi « 27r x 1 Hz, is also listed 
for comparison, along with our anticipated experimental sen- 
sitivity (corresponding to 5il = 27r x 10 fjJiz in 1 second of 
integration). An integration time of 10® seconds is assumed. 
For fc = 1, Xp = fi/c, leading to a spin-gravity signal at the 
level AK«3 x 10" ^^ 



Description 



A3? 



Atomic shot-noise limit 4 x 10^^^ 

Anticipated sensitivity 1 x 10~^^ 

Gyro-compass effect 2 x 10^^^ 

Nuclear magnetic moments < 10~^® 

Nonlinear Zeeman effect negligible 

Light shifts < 10"^^ 

Spin-exchange collisions 2 x 10^^* 

Magnetic field gradients & geometric phase < 10~^® 

Wall collisions < 10'^*^ 

effect, known as the gyro-compass effect [33] or the spin- 
rotation effect [ij, can be understood as the result of 
viewing an inertial system, the atomic spins, from a non- 
inertial frame, the surface of the rotating Earth. Uncer- 
tainty in the magnitude of this i?-odd systematic effect 
can be made quadratic in the misalignment of the ex- 
perimental apparatus (Eq. [25|) by orienting B along the 
axis of Earth's rotation CIe [311]. This approach has been 
implemented as shown in the picture at the top left of 
the experimental setup diagram (Fig. [5]) . 

Including the gyro-compass effect adds a _B-odd spin- 
precession frequency {Clo = flECOsB) to fls5 and J^gy, 
where 6 describes the misalignment between Earth's ro- 
tation axis fls and B. Based on Eq. p6|) : 



A3i 



787 - 785 
787 + 785 



4.12 



Xp9 cos ( 



-2.4 



He cos 9 



(26) 



The angle is now the resultant angle between g and 
flE {(j) equals 90° plus the latitude of the laboratory lo- 
cation, about 37°, so coscj) ~ —0.6). We can control the 
orientation of B with respect to an auxiliary laser beam 
propagating along z to a level of better than 10~^ [5a |. 
The long lever arm of the laser beam (in combination 
with GPS and aerial surveying) enables alignment of the 
auxiliary laser beam propagation direction with Qe to 
within « 0.3° « 5 x 10~'^ rad, so that systematic uncer- 
tainty in the gyro-compass effect due apparatus misalign- 
ment is at the 3 x 10~^° Hz level. Thus errors due to the 
Earth's rotation can be well-controlled at our proposed 
level of sensitivity. 



C. Nuclear magnetic moments 

Although nuclear magnetic moments are a thousand 
times smaller than i^q, their effect on the observed spin 
precession frequencies is clearly evident in the sample 
data shown in Fig. |3l giving rise to a slow beating visible 
in the time-domain signal shown in the upper plot. The 
nuclear magnetic moment modifies the Lande factors for 
the alkali ground state hyperfine levels [6l, |6a] : 



9F=I+i 



9f=i- 



2/+1 



91 



fJ-N 



21 



^io 21 + 1 



2/+1 



91- 



1) 



Mo 



2/-fl 



(27) 
(28) 



where gj is the nuclear Lande factor {gj sa 0.539 for *^Rb, 
gj w 1.827 for *'^Rb [63) Mat is the nuclear magneton, 
and /^at/mo ~ 5 x 10~^. This creates a difference in the 
Larmor frequencies for atoms in the two different ground 
state hyperfine levels 



= -2gifJ.NB . 



l]-n. 



2 



(29) 



For B w 7.1443 mG as in the data shown in Fig.|31 *^Rb 
hasAfinuci ~ — 27rx 5.9 Hz and ^"^ Kb has Annuel ~ — 27rx 
19.9 Hz. Off-resonant synchronous optical pumping for 
the *^Rb F — 2 state is much more efficient than that 
for the *^Rb F = 1 state because of the smaller Ailnuci. 
The signal from the ^^Rb F ^ 1 state is additionally 
suppressed relative to the signal from the ^'^Rb F — 2 
state because the re-pump laser beam depletes the *^Rb 
F — 1 state and the probe laser light is farther detuned 
from the ^^Rb D2 F = 1 —>■ F' resonance than from the 
*^Rb D2 F = 2 —i' F' resonance. Consequently, only the 
dominant signals from the '^^Rb F — 3 and *^Rb F — 2 
states, along with a much smaller signal from the *^Rb 
F = 2 state which leads to the slow beating observed in 
the time-domain signal in Fig. [3l are easily detectable in 
the data. 

In our data analysis, the three observable resonances 
in the Fourier transformed optical rotation data are fit 
directly, and modeling demonstrates that neglecting the 
resonance associated with the ^^Rb F = I state in our 
fitting routine does not affect our analysis at the desired 
level of accuracy. Fortunately, any first-order system- 
atic effect associated with the nuclear magnetic moments 
manifests as a i?-even systematic that is suppressed by 
~ 5B/B ~ 10~^ in the comagnetometer signal A3J as 
discussed in Sec. IV Al 



D. Nonlinear Zeeman effect 

The magnetic field also mixes Zeeman sublevels in dif- 
ferent ground state hyperfine levels, leading the Zeeman 
effect to acquire a nonlinear dependence on B. In our 
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experiment, the nonlinear Zeeman effect manifests as a 
splitting of the Larmor resonances [68. - . 71.] . The splitting 
of the resonances is symmetric about the unperturbed 
Larmor frequency and smaller than the linewidth, there- 
fore, to leading order, it does not contribute any system- 
atic shift to the spin precession frequencies. The energy 
E{F,Mp) of a particular ground state Zeeman sublevel 
{Mp is the projection of F along z) of an alkali atom is 
described by the Breit-Rabi formula |72| : 

E{F ^ I ±1/2, Mp) = 



A 



hfs 



A 



hfs , 



^ g,^^BMp±^-^{2I+l)Jl + ^ 



(30) 



where ^Ihfs is the alkali atom's hyperfine structure con- 
stant and u is the perturbation parameter given by: 



gj^-o + ff/MJV 2B 



HqB 



2/+1 A 



hfs 



2/ + iyi 



hfs 



(31) 



where gj « 2 is the Lande g-factor for the electron. The 
Breit-Rabi formula (Eq. [50)1 can be expanded to second 
order in u and the terms proportional to u'^ can be iden- 
tified as the nonlinear Zeeman shifts E^W- 



EnUF^I± 1/2, Mp 

2-^hfs 



±u' 



(2/+1) 1- 



4Af| 



(2/ +1)2 
The term in Eq. (|32|) proportional to Mp, 



(32) 



T 



^^lB^ 



{21 + ly ^hfs 



M 



F ' 



(33) 



causes a nonlinear Zeeman shift of the Larmor frequen- 
cies that splits a single spin-precession resonance into 
multiple resonances. In the following we consider only 
the F = 1 + 1/2 ground state hyperfine levels and de- 
fine the unperturbed Larmor frequency ft)^ as the term 
linear in B, 

^i°^ - (^J-^fio-9if^N:^^)B . (34) 



2/ -1-1' 



21 + 1 



For the *^Rb F = 3 state, there appear six resonance 



frequencies split symmetrically about flj^ 

(21+ if ^hfs ' 

(21+ if Ams ' 

8 M§i3^ 
(2/+l)^ ^hfs 



(0). 



fif ±5 



fif ± 3 



l]f ± 1 



(For the *^Rb F = 2 state there are four resonance fre- 
quencies described by the latter four cases above.) For 



B sa 7.1443 mG, as in the data shown in Fig.|3l the max- 
imum splitting of the resonance frequencies is ~ 0.2 Hz 
for the ^5Rb F ^ 3 state and w 0.08 Hz for the ^'^Rb 
F = 2 state, in both cases smaller than the resonance 
linewidth of ~ 1 Hz. Any imbalance in the population 
of the Zeeman sublevels associated with the different res- 
onances constitutes longitudinal spin polarization that 
does not contribute to the spin precession signal. There- 
fore, in some sense, the signal amplitudes for different 
resonances are naturally balanced. Thus the only ap- 
parent consequence of the nonlinear Zeeman effect under 
our experimental conditions is a slight broadening of the 
spin precession resonances. Nonetheless, measurements 
will be carried out at different magnetic fields to test for 
any magnetic-field-dependent systematic errors. 



E. Light shifts 

The ac Stark effect due to the optical electric field of 
the probe beam can cause light shifts of Zeeman sub- 
levels, leading to shifts of the measured precession fre- 
quencies for ^^Rb and ®^Rb. In general, ac Stark shifts 
can be described in terms of scalar, vector, and tensor 
polarizabilities [6a, [Tj]. The scalar and tensor polariz- 
abilities are described by rank-zero and rank-two opera- 
tors, and thus their effect on atoms can be modeled as a 
fictitious static electric field along the light polarization 
axis; the vector polarizability is described by a rank-one 
operator, and thus can be modeled as a fictitious static 
ma gne tic field along the light propagation direction k 

0m- 

For measurement of ilss a-nd rjg?, it is the vector light 
shift in particular that causes the most significant sys- 
tematic effect. Although the probe beam is nominally 
linearly polarized and detuned far from the Doppler- 
broadened optical resonances, vector light shifts can still 
arise due to residual ellipticity e induced in the beam due 
to birefringence of the vapor cell walls. Measurements of 
the probe beam polarization before and after the cell us- 
ing a ThorLABs PAX720IR1-T polarimeter system show 



that e can be made < 0.01° « 2 x 10" 



rad. 



The frequency shift fiac associated with the vector po- 
larizability is a i?-odd systematic effect, and thus is not 
suppressed in the comagnetometer signal A!R. However, 
because k is orthogonal to B, there is a geometric sup- 
pression according to Eq. ([^^ . The quadratic correction 
term in Eq. ((22)) can be neglected in our case, and we 
have: 



J7ac ~ v9sin(2e)ASac/ft ~ 2ipeAEi^c/fi 



(35) 



where A_Bac is the vector light shift between adjacent 
Zeeman sublevels {AMp = 1) for left-circularly polarized 
light along the quantization axis. AE^c can be estimated, 
for example, based on the formula from Ref. 



M- 



AE, 



|(5^i/2|er|5Pi/2) 



9Aw, 



3/2 



-3i=^(|£o| 



(36) 
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where (5S'i/2|er|5Pi/2) ~ 3eao is the transition dipole 
matrix element between the bSi/2 and 5Pi /2 states, Oq is 
the Bohr radius, AcL'3/2 is the detuning of the probe beam 
from the D2 resonance, gp is the ground state Lande 
factor, and (|£oP) is the average square of the optical 
electric field experienced by the atoms. In calculating 
(|£oP), one must take into account the fact that the ef- 
fective optical electric field experienced by the atoms is 
diluted by the ratio of the volume within the cell illu- 
minated by the probe light beam to the total volume 
of the cell (for our experiment, the ratio « 2 x 10"'^) 
since the atoms spend only a small fraction of their time 
in the probe light during the precession time [70|- For 
our typical probe light power of 200 /xW and detuning 
of w 3 GHz below the center frequency of the Doppler- 
broadened ^^Rb D2 F = 2 —i' F' resonance, we estimate 
that the vector light shifts for ^'"'Rb and *^Rb are given, 
respectively, by 



sion rate 7so <^ f^L, we have Ref. |80| : 



n. 



Is 



AE^c{S5) « -2Trh x (3 Hz) , 
A£;ac(87) « -2TTh X (1.3 Hz) 



(37) 
(38) 



Because sensitive nonlinear magneto-optical effects can 
be used to directly measure the angle between k and 
B [56|, it is actually feasible in our setup to constrain 
if < 10"^ rad. Therefore, based on Eq. (I3S|), fiac < 
10~® Hz, and consequently light shifts are not expected 
to prevent the experiment from reaching its sensitivity 
target. Nevertheless, data will be taken at different probe 
light powers to check for any systematic effects related to 
light shifts. 



F. Spin-exchange collisions 

The dual-isotope Rb comagnetometer relies on inde- 
pendent measurements of Slgs and flsTj so coupling be- 
tween the two isotopes through spin-exchange (SE) colli- 
sions can produce a systematic error. However, since the 
experiment is carried out in a bias field of |B| ~ 10 mG 
and rigs 7^ i^87: in the frame rotating with each isotope's 
precession frequency, the spin-polarization of the other 
isotope is time-averaged to nearly zero Q- Nonethe- 
less, there still appears a small SE frequency shift [80| . 
Spin-exchange collisions tend to pull the precession fre- 
quencies toward a weighted average: SE collisions that 
transfer atoms between ground state hyperfine levels of 
a single isotope reduce the measured spin precession fre- 
quency since the gyromagnetic ratios have opposite signs; 
SE collisions between Rb isotopes shift figs to a higher 
frequency and flsr to a lower frequency (cross-isotope 
SE shifts do not cancel because of the larger statistical 
weights of the F = I + 1/2 hyperfine levels). 

An estimate of the scale of the SE frequency shift ilgo 
can be obtained by considering SE collisions between 
ground state hyperfine levels of each individual isotope. 
Under our experimental conditions, where the SE colli- 



1 



1 



1 



IS^JlV^ (2/+1)W^ (2/+l)V ■ ^^^^ 
Under the experimental conditions for the data shown 
in Fig. 13 (7sc ^ 2n x 1.3 Hz, flgz w 27r x 3334 Hz, and 
risT w 27r X 5001 Hz): 

f^so(85) « -2tt X 2.3 X 10"^ Hz , (40) 

flsciSl) w ~2n X 1.3 X 10"^ Hz . (41) 

Crucially, SE frequency shifts are i?-even and so their 
effect on the comagnetometer signal A3? is described by 
Eq. (|17l) . thus suppressing any SE collision-related sys- 
tematic effects by ^ SB/ B ^ 10^^. Because of the sup- 
pression of Vise in A3?, systematic effects due to SE colli- 
sions between atoms are negligible in our experiment. 



G. Other systematic effects 

Another concern is the effect of magnetic field gradi- 
ents which cause *^Rb and *^Rb atoms to, on a ran- 
dom basis, sample different magnetic fields, reducing 
the effectiveness of the comagnetometry scheme. Field 
gradients are nulled using auxiliary measurements to 
< 10^^ G/cm in all directions |50j. Effects of gradi- 
ents are further reduced due to motional averaging in 
the evacuated antirelaxation-coated cells: atoms typi- 
cally bounce off of the cell walls ^10^ times between 
interactions with the laser beam [44]. For a sample of 
^ 10^^ atoms, this creates uncertainty at the uHz/a/Hz 
level, well below our statistical sensitivity to spin preces- 
sion. Furthermore, systematic frequency shifts related 
to the geometric (Berry's) phase [81| are proportional to 
gradients, and are estimated to be less than a nHz under 
typical experimental conditions based on the analysis of 
Ref. [13. 

Wall collisions can produce quadrupolar splittings of 
spin precession frequencies due interaction of atomic 
spins with surface electric field gradients [3l|. In our ex- 
periment, wall collisions should produce negligibly small 
shifts of ilss E^nd Jig? since the vapor cell employs an 
amorphous antirelaxation coating and is spherical in 
shape, so that, to a high precision, there is no preferred 
direction in the cell. We can estimate that in the worst- 
case scenario the contribution to a cell-related shift is 
on the order of the wall relaxation rate (~ 10~^ Hz) 
times the square of the ratio of the size of the open- 
ing to the stem that contains the alkali metal sam- 
ple (« 10~^ cm^) to the inner surface area of the cell 
(« 80 cm2): < 10^1° Hz. 



VI. CONCLUSION 

An experiment measuring spin precession frequencies 
of overlapping ensembles of ^^Rb and *^Rb atoms con- 
tained within an evacuated, antirelaxation-coated vapor 
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cell can be used to search for presently unconstrained 
anomalous long-range spin-mass couplings. Synchronous 
optical pumping with circularly polarized light is used 
to generate spin polarization transverse to an applied 
magnetic field and optical rotation of a linearly polarized 
probe beam is used to measure the ®^Rb and ^^Rb spin 
precession frequencies. The Earth is used as the source 
mass. The present statistical sensitivity of the apparatus 
to spin precession frequencies is 10~^ Hz in one second 
of integration, with a shot-noise-projected sensitivity ex- 
ceeding this level by over two orders of magnitude. A va- 
riety of systematic errors are considered, and all known 
sources of error can be controlled at the 10~^ Hz level. 

There are several promising and potentially more sen- 
sitive approaches to searching for long-range spin-mass 
couplings, including the use of spin-exchange- relaxation 
free (SERF) comagnetometers (sljp [H, ^Be/^^^Xe 
free-precession comagnetometers [831 , and liquid state 
nuclear-spin comagnetometers [86j. However, there are 
experimental challenges to applying each of these alter- 
native approaches to a search for long-range spin-mass 
couplings. For example, it is potentially difficult to dis- 
tinguish the coupling of spins to the local gravitational 
field from other lab-fixed backgrounds with a SERF co- 



magnetometer [83|, and magnetic field gradients may be 
an issue for liquid state nuclear-spin comagnetometers 

These same techniques can also be applied to search for 
long-range anomalous spin-spin interactions using polar- 
ized electrons in the Earth [83. 
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